Duet: SNP-assisted structural variant calling and phasing using Oxford nanopore sequencing.

BACKGROUND: Whole genome sequencing using the long-read Oxford Nanopore Technologies (ONT) MinION sequencer provides a cost-effective option for structural variant (SV) detection in clinical applications. Despite the advantage of using long reads, however, accurate SV calling and phasing are still challenging. RESULTS: We introduce Duet, an SV detection tool optimized for SV calling and phasing using ONT data. The tool uses novel features integrated from both SV signatures and single-nucleotide polymorphism signatures, which can accurately distinguish SV haplotype from a false signal. Duet was benchmarked against state-of-the-art tools on multiple ONT sequencing datasets of sequencing coverage ranging from 8× to 40×. At low sequencing coverage of 8×, Duet performs better than all other tools in SV calling, SV genotyping and SV phasing. When the sequencing coverage is higher (20× to 40×), the F1-score for SV phasing is further improved in comparison to the performance of other tools, while its performance of SV genotyping and SV calling remains higher than other tools. CONCLUSION: Duet can perform accurate SV calling, SV genotyping and SV phasing using low-coverage ONT data, making it very useful for low-coverage genomes. It has great performance when scaled to high-coverage genomes, which is adaptable to various clinical applications. Duet is open source and is available at https://github.com/yekaizhou/duet .

Clair3-trio: high-performance Nanopore long-read variant calling in family trios with trio-to-trio deep neural networks.

Accurate identification of genetic variants from family child-mother-father trio sequencing data is important in genomics. However, state-of-the-art approaches treat variant calling from trios as three independent tasks, which limits their calling accuracy for Nanopore long-read sequencing data. For better trio variant calling, we introduce Clair3-Trio, the first variant caller tailored for family trio data from Nanopore long-reads. Clair3-Trio employs a Trio-to-Trio deep neural network model, which allows it to input the trio sequencing information and output all of the trio's predicted variants within a single model to improve variant calling. We also present MCVLoss, a novel loss function tailor-made for variant calling in trios, leveraging the explicit encoding of the Mendelian inheritance. Clair3-Trio showed comprehensive improvement in experiments. It predicted far fewer Mendelian inheritance violation variations than current state-of-the-art methods. We also demonstrated that our Trio-to-Trio model is more accurate than competing architectures. Clair3-Trio is accessible as a free, open-source project at https://github.com/HKU-BAL/Clair3-Trio.

Identification of a Novel Betacoronavirus (Merbecovirus) in Amur Hedgehogs from China.

While dromedaries are the immediate animal source of Middle East Respiratory Syndrome (MERS) epidemic, viruses related to MERS coronavirus (MERS-CoV) have also been found in bats as well as hedgehogs. To elucidate the evolution of MERS-CoV-related viruses and their interspecies transmission pathway, samples were collected from different mammals in China. A novel coronavirus related to MERS-CoV, Erinaceus amurensis hedgehog coronavirus HKU31 (Ea-HedCoV HKU31), was identified from two Amur hedgehogs. Genome analysis supported that Ea-HedCoV HKU31 represents a novel species under Merbecovirus, being most closely related to Erinaceus CoV from European hedgehogs in Germany, with 79.6% genome sequence identity. Compared to other members of Merbecovirus, Ea-HedCoV HKU31 possessed unique non-structural proteins and putative cleavage sites at ORF1ab. Phylogenetic analysis showed that Ea-HedCoV HKU31 and BetaCoV Erinaceus/VMC/DEU/2012 were closely related to NeoCoV and BatCoV PREDICT from African bats in the spike region, suggesting that the latter bat viruses have arisen from recombination between CoVs from hedgehogs and bats. The predicted HKU31 receptor-binding domain (RBD) possessed only one out of 12 critical amino acid residues for binding to human dipeptidyl peptidase 4 (hDPP4), the MERS-CoV receptor. The structural modeling of the HKU31-RBD-hDPP4 binding interphase compared to that of MERS-CoV and Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) suggested that HKU31-RBD is unlikely to bind to hDPP4. Our findings support that hedgehogs are an important reservoir of Merbecovirus, with evidence of recombination with viruses from bats. Further investigations in bats, hedgehogs and related animals are warranted to understand the evolution of MERS-CoV-related viruses.

Replication of MERS and SARS coronaviruses in bat cells offers insights to their ancestral origins.

Previous findings of Middle East Respiratory Syndrome coronavirus (MERS-CoV)-related viruses in bats, and the ability of Tylonycteris-BatCoV HKU4 spike protein to utilize MERS-CoV receptor, human dipeptidyl peptidase 4 hDPP4, suggest a bat ancestral origin of MERS-CoV. We developed 12 primary bat cell lines from seven bat species, including Tylonycteris pachypus, Pipistrellus abramus and Rhinolophus sinicus (hosts of Tylonycteris-BatCoV HKU4, Pipistrellus-BatCoV HKU5, and SARS-related-CoV respectively), and tested their susceptibilities to MERS-CoVs, SARS-CoV, and human coronavirus 229E (HCoV-229E). Five cell lines, including P. abramus and R. sinicus but not T. pachypus cells, were susceptible to human MERS-CoV EMC/2012. However, three tested camel MERS-CoV strains showed different infectivities, with only two strains capable of infecting three and one cell lines respectively. SARS-CoV can only replicate in R. sinicus cells, while HCoV-229E cannot replicate in any bat cells. Bat dipeptidyl peptidase 4 (DPP4) sequences were closely related to those of human and non-human primates but distinct from dromedary DPP4 sequence. Critical residues for binding to MERS-CoV spike protein were mostly conserved in bat DPP4. DPP4 was expressed in the five bat cells susceptible to MERS-CoV, with significantly higher mRNA expression levels than those in non-susceptible cells (P = 0.0174), supporting that DPP4 expression is critical for MERS-CoV infection in bats. However, overexpression of T. pachypus DPP4 failed to confer MERS-CoV susceptibility in T. pachypus cells, suggesting other cellular factors in determining viral replication. The broad cellular tropism of MERS-CoV should prompt further exploration of host diversity of related viruses to identify its ancestral origin.

Discovery and Sequence Analysis of Four Deltacoronaviruses from Birds in the Middle East Reveal Interspecies Jumping with Recombination as a Potential Mechanism for Avian-to-Avian and Avian-to-Mammalian Transmission.

The emergence of Middle East respiratory syndrome showed once again that coronaviruses (CoVs) in animals are potential source for epidemics in humans. To explore the diversity of deltacoronaviruses in animals in the Middle East, we tested fecal samples from 1,356 mammals and birds in Dubai, The United Arab Emirates. Four novel deltacoronaviruses were detected from eight birds of four species by reverse transcription-PCR (RT-PCR): FalCoV UAE-HKU27 from a falcon, HouCoV UAE-HKU28 from a houbara bustard, PiCoV UAE-HKU29 from a pigeon, and QuaCoV UAE-HKU30 from five quails. Complete genome sequencing showed that FalCoV UAE-HKU27, HouCoV UAE-HKU28, and PiCoV UAE-HKU29 belong to the same CoV species, suggesting recent interspecies transmission between falcons and their prey, houbara bustards and pigeons, possibly along the food chain. Western blotting detected specific anti-FalCoV UAE-HKU27 antibodies in 33 (75%) of 44 falcon serum samples, supporting genuine infection in falcons after virus acquisition. QuaCoV UAE-HKU30 belongs to the same CoV species as porcine coronavirus HKU15 (PorCoV HKU15) and sparrow coronavirus HKU17 (SpCoV HKU17), discovered previously from swine and tree sparrows, respectively, supporting avian-to-swine transmission. Recombination involving the spike protein is common among deltacoronaviruses, which may facilitate cross-species transmission. FalCoV UAE-HKU27, HouCoV UAE-HKU28, and PiCoV UAE-HKU29 originated from recombination between white-eye coronavirus HKU16 (WECoV HKU16) and magpie robin coronavirus HKU18 (MRCoV HKU18), QuaCoV UAE-HKU30 from recombination between PorCoV HKU15/SpCoV HKU17 and munia coronavirus HKU13 (MunCoV HKU13), and PorCoV HKU15 from recombination between SpCoV HKU17 and bulbul coronavirus HKU11 (BuCoV HKU11). Birds in the Middle East are hosts for diverse deltacoronaviruses with potential for interspecies transmission.IMPORTANCE During an attempt to explore the diversity of deltacoronaviruses among mammals and birds in Dubai, four novel deltacoronaviruses were detected in fecal samples from eight birds of four different species: FalCoV UAE-HKU27 from a falcon, HouCoV UAE-HKU28 from a houbara bustard, PiCoV UAE-HKU29 from a pigeon, and QuaCoV UAE-HKU30 from five quails. Genome analysis revealed evidence of recent interspecies transmission between falcons and their prey, houbara bustards and pigeons, possibly along the food chain, as well as avian-to-swine transmission. Recombination, which is known to occur frequently in some coronaviruses, was also common among these deltacoronaviruses and occurred predominantly at the spike region. Such recombination, involving the receptor binding protein, may contribute to the emergence of new viruses capable of infecting new hosts. Birds in the Middle East are hosts for diverse deltacoronaviruses with potential for interspecies transmission.

Receptor Usage of a Novel Bat Lineage C Betacoronavirus Reveals Evolution of Middle East Respiratory Syndrome-Related Coronavirus Spike Proteins for Human Dipeptidyl Peptidase 4 Binding.

Although bats are known to harbor Middle East Respiratory Syndrome coronavirus (MERS-CoV)-related viruses, the role of bats in the evolutionary origin and pathway remains obscure. We identified a novel MERS-CoV-related betacoronavirus, Hp-BatCoV HKU25, from Chinese pipistrelle bats. Although it is closely related to MERS-CoV in most genome regions, its spike protein occupies a phylogenetic position between that of Ty-BatCoV HKU4 and Pi-BatCoV HKU5. Because Ty-BatCoV HKU4 but not Pi-BatCoV HKU5 can use the MERS-CoV receptor human dipeptidyl peptidase 4 (hDPP4) for cell entry, we tested the ability of Hp-BatCoV HKU25 to bind and use hDPP4. The HKU25-receptor binding domain (RBD) can bind to hDPP4 protein and hDPP4-expressing cells, but it does so with lower efficiency than that of MERS-RBD. Pseudovirus assays showed that HKU25-spike can use hDPP4 for entry to hDPP4-expressing cells, although with lower efficiency than that of MERS-spike and HKU4-spike. Our findings support a bat origin of MERS-CoV and suggest that bat CoV spike proteins may have evolved in a stepwise manner for binding to hDPP4.

Bats host diverse parvoviruses as possible origin of mammalian dependoparvoviruses and source for bat-swine interspecies transmission.

Compared to the enormous species diversity of bats, relatively few parvoviruses have been reported. We detected diverse and potentially novel parvoviruses from bats in Hong Kong and mainland China. Parvoviruses belonging to Amdoparvovirus, Bocaparvovirus and Dependoparvovirus were detected in alimentary, liver and spleen samples from 16 different chiropteran species of five families by PCR. Phylogenetic analysis of partial helicase sequences showed that they potentially belonged to 25 bocaparvovirus, three dependoparvovirus and one amdoparvovirus species. Nearly complete genome sequencing confirmed the existence of at least four novel bat bocaparvovirus species (Rp-BtBoV1 and Rp-BtBoV2 from Rhinolophus pusillus, Rs-BtBoV2 from Rhinolophus sinicus and Rol-BtBoV1 from Rousettus leschenaultii) and two novel bat dependoparvovirus species (Rp-BtAAV1 from Rhinolophus pusillus and Rs-BtAAV1 from Rhinolophus sinicus). Rs-BtBoV2 was closely related to Ungulate bocaparvovirus 5 with 93, 72.1 and 78.7 % amino acid identities in the NS1, NP1 and VP1/VP2 genes, respectively. The detection of bat bocaparvoviruses, including Rs-BtBoV2, closely related to porcine bocaparvoviruses, suggests recent interspecies transmission of bocaparvoviruses between bats and swine. Moreover, Rp-BtAAV1 and Rs-BtAAV1 were most closely related to human AAV1 with 48.7 and 57.5 % amino acid identities in the rep gene. The phylogenetic relationship between BtAAVs and other mammalian AAVs suggests bats as the ancestral origin of mammalian AAVs. Furthermore, parvoviruses of the same species were detected from multiple bat species or families, supporting the ability of bat parvoviruses to cross species barriers. The results extend our knowledge on the diversity of bat parvoviruses and the role of bats in parvovirus evolution and emergence in humans and animals.

Identification and interspecies transmission of a novel bocaparvovirus among different bat species in China.

We report the discovery of a novel bocaparvovirus, bat bocaparvovirus (BtBoV), in one spleen, four respiratory and 61 alimentary samples from bats of six different species belonging to three families, Hipposideridae, Rhinolophidae and Vespertilionidae. BtBoV showed a higher detection rate in alimentary samples of Rhinolophus sinicus (5.7 %) than those of other bat species (0.43-1.59 %), supporting R. sinicus as the primary reservoir and virus spillover to accidental bat species. BtBoV peaked during the lactating season of R. sinicus, and it was more frequently detected among female than male adult bats (P<0.05), and among lactating than non-lactating female bats (P<0.0001). Positive BtBoV detection was associated with lower body weight in lactating bats (P<0.05). Ten nearly complete BtBoV genomes from three bat species revealed a unique large ORF1 spanning NS1 and NP1 in eight genomes and conserved splicing signals leading to multiple proteins, as well as a unique substitution in the conserved replication initiator motif within NS1. BtBoV was phylogenetically distantly related to known bocaparvoviruses with ≤57.3 % genome identities, supporting BtBoV as a novel species. Ms-BtBoV from Miniopterus schreibersii and Hp-BtBoV from Hipposideros pomona demonstrated 97.2-99.9 % genome identities with Rs-BtBoVs from R. sinicus, supporting infection of different bat species by a single BtBoV species. Rs-BtBoV_str15 represents the first bat parvovirus genome with non-coding regions sequenced, which suggested the presence of head-to-tail genomic concatamers or episomal forms of the genome. This study represents the first to describe interspecies transmission in BoVs. The high detection rates in lactating female and juvenile bats suggest possible vertical transmission of BtBoV.

Severe Acute Respiratory Syndrome (SARS) Coronavirus ORF8 Protein Is Acquired from SARS-Related Coronavirus from Greater Horseshoe Bats through Recombination.

UNLABELLED: Despite the identification of horseshoe bats as the reservoir of severe acute respiratory syndrome (SARS)-related coronaviruses (SARSr-CoVs), the origin of SARS-CoV ORF8, which contains the 29-nucleotide signature deletion among human strains, remains obscure. Although two SARS-related Rhinolophus sinicus bat CoVs (SARSr-Rs-BatCoVs) previously detected in Chinese horseshoe bats (Rhinolophus sinicus) in Yunnan, RsSHC014 and Rs3367, possessed 95% genome identities to human and civet SARSr-CoVs, their ORF8 protein exhibited only 32.2 to 33% amino acid identities to that of human/civet SARSr-CoVs. To elucidate the origin of SARS-CoV ORF8, we sampled 348 bats of various species in Yunnan, among which diverse alphacoronaviruses and betacoronaviruses, including potentially novel CoVs, were identified, with some showing potential interspecies transmission. The genomes of two betacoronaviruses, SARSr-Rf-BatCoV YNLF_31C and YNLF_34C, from greater horseshoe bats (Rhinolophus ferrumequinum), possessed 93% nucleotide identities to human/civet SARSr-CoV genomes. Although these two betacoronaviruses displayed lower similarities than SARSr-Rs-BatCoV RsSHC014 and Rs3367 in S protein to civet SARSr-CoVs, their ORF8 proteins demonstrated exceptionally high (80.4 to 81.3%) amino acid identities to that of human/civet SARSr-CoVs, compared to SARSr-BatCoVs from other horseshoe bats (23.2 to 37.3%). Potential recombination events were identified around ORF8 between SARSr-Rf-BatCoVs and SARSr-Rs-BatCoVs, leading to the generation of civet SARSr-CoVs. The expression of ORF8 subgenomic mRNA suggested that the ORF8 protein may be functional in SARSr-Rf-BatCoVs. The high Ka/Ks ratio among human SARS-CoVs compared to that among SARSr-BatCoVs supported that ORF8 is under strong positive selection during animal-to-human transmission. Molecular clock analysis using ORF1ab showed that SARSr-Rf-BatCoV YNLF_31C and YNLF_34C diverged from civet/human SARSr-CoVs in approximately 1990. SARS-CoV ORF8 originated from SARSr-CoVs of greater horseshoe bats through recombination, which may be important for animal-to-human transmission. IMPORTANCE: Although horseshoe bats are the primary reservoir of SARS-related coronaviruses (SARSr-CoVs), it is still unclear how these bat viruses have evolved to cross the species barrier to infect civets and humans. Most human SARS-CoV epidemic strains contain a signature 29-nucleotide deletion in ORF8, compared to civet SARSr-CoVs, suggesting that ORF8 may be important for interspecies transmission. However, the origin of SARS-CoV ORF8 remains obscure. In particular, SARSr-Rs-BatCoVs from Chinese horseshoe bats (Rhinolophus sinicus) exhibited <40% amino acid identities to human/civet SARS-CoV in the ORF8 protein. We detected diverse alphacoronaviruses and betacoronaviruses among various bat species in Yunnan, China, including two SARSr-Rf-BatCoVs from greater horseshoe bats that possessed ORF8 proteins with exceptionally high amino acid identities to that of human/civet SARSr-CoVs. We demonstrated recombination events around ORF8 between SARSr-Rf-BatCoVs and SARSr-Rs-BatCoVs, leading to the generation of civet SARSr-CoVs. Our findings offer insight into the evolutionary origin of SARS-CoV ORF8 protein, which was likely acquired from SARSr-CoVs of greater horseshoe bats through recombination.